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Abstract 

The Ostrogradsky instability of higher derivative Lagrangians is de- 
rived from first principles using Control theory and Lyapunov Stability 
Analysis. This result is then used to argue that Born-Infeld Lagrangians 
are viable modifications of the Einstein-Hilbert action provided the action 
is varied in accordance with the Palatini formalism, in contrast to the met- 
ric formalism. Finally, the Born-Infeld version of the FRW equations are 
derived and the cosmological dynamics is studied for matter dominated 
closed and open universes, and the results are compared with the usual 
cosmology. 

Introduction 

The light curves plotted from several hundred Type la supernovae [1] indicate 
that the universe is expanding at an accelerated rate. On the other hand, the 
Friedmann equations derived from Einstein's General Relativity show that for a 
fluid that sat- isfies the Weak Energy Condition, the universe should decelerate 
as it expands [2]. This astonishing discrepancy between observation and the 
prediction of the general theory of relativity led to a spurt of research activities 
intended to modify the existing cosmological models with an aim to make them 
consistent with observation. 

The endeavor started with models that include either the Cosmological Con- 
stant, Quintessence, Dark Energy, or Phantom Energy [3] that presumably 
drives the accelerated expansion. However, all such models faced one or the 
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other problems. For instance, in the case of the models with cosmological con- 
stant, the smallness of its value inferred from astronomical data compared to 
its value predicted by particle physics and the remarkable but implausible co- 
incidence of the vacuum energy density and the mass density in the present 
epoch, require two-fold fine tuning which does not appear appealing at all. On 
the other hand, models advocating for the existence of dark energy suffer from 
their own problems related to the very nature of dark energy and the fact that 
different models predict different rates of expansion that require a rather tedious 
task of searching for the fingerprints of the dark energy in the entire history of 
cosmic expansion. Added to it is the unpleasant fact that distinguishing among 
different histories of expansion would require measurements an order of magni- 
tude more accurate than the existing ones [12]. 

Apart from that, in order to make Einstein's field equation with dark energy 
agree with the observed anisotropics in the cosmic microwave background, the 
amount of dark matter that should be posited is eighteen times the observed 
ordinary matter. The unresolved, so far, problems related to these models 
inspired a third route, namely, the modification of the Einstein's field equa- 
tions themselves by taking into account the higher order derivative terms in the 
Einstein-Hilbert (E-H) action. 

However, as it turned out, this cannot be accomplished in an offhand man- 
ner. There is a no-go theorem that constrains the form of the Lagrangian that 
can be used for this purpose. In its essence, it is related to Newton's observa- 
tion that the equations of physics, when written in terms of the fundamental 
quantities, are always second-order in time [3]. This no-go theorem is called Os- 
trogradsky's theorem: "A system whose Lagrangian depends non-degenerately 
on the second- and higher-order derivatives of the dynamical quantities is nec- 
essarily un- stable". The terminology non-degenerately means that the higher 
derivatives can be solved in terms of the lower derivatives. This theorem has 
been invoked in the literature to support the claim that the modifications of the 
E-H Lagrangian that depend on the traces of the Riemann and Ricci tensors 
must be excluded from consideration [3]. But the justification of this important 
result has been done by bringing ideas of second quantization into the classical 
setting. Moreover, the sign indefiniteness of the energy only shows that the 
energy function is not a Lyapunov function. It does not indicate that no other 
Lyapunov function can exist. 

In the given paper, we first show that the proof of Ostrogradsky's theorem 
can be made more precise within classical mechanics, without discussing the 
particle picture. Then we proceed on to show that, if we consider a f(R) theory 
of gravity by modifying the E-H Lagrangian with the Born-Infeld (B-I) type of 
nonlinearity, Ostrogradsky's theorem can be circumvented under the condition 
that the modified action is varied only in Palatini's formalism. In view of this, 
our work is devoted to the study of the B-I modified gravity model in Palatini's 
formalism. The cosmological implications of such a model are studied for the 
case when the spatial hyper-surface can be curved. Note that most of the works 
on B-I cosmology [11] have been restricted to the case of flat spatial hypersur- 
faces. 
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The organization of the paper is as follows. We start with the proof of 
Ostrogradsky's theorem in the classical mechanical framework. We then take 
up the modified E-H Lagrangian with B-I nonlinearity and justify why the 
corresponding action should be varied only in the Palatini's formalism. Finally, 
we derive Friedmann-Robertson- Walker (FRW) equations for the B-I gravity, 
study their cosmological implications and present the results of our numerical 
calculations 

Chaetev's Theorem 

We first review a theorem from control theory. This result is very important 
for the subsequent discussion with regard to the proof of the Ostrogradsky 
instability. 

Consider the following system of differential equations 

x = /(x) (1) 

Here, x represents an n-dimensional vector, while f(x) is a scalar point function. 
Without losing generality, we can assume that x = is the equilibrium point of 
the system. We now state the result. The proof can be looked up in [B]. 

Chaetev's Theorem [6] : Let x = be an equilibrium point of the system 
for (1). Let V : D — > R be a continously differentiable function such that 
V(0) = and ^(xq) > for some x with ||x || arbitrarily small. For r > 0, 
let B r denote the set of all x with | |x| | less than r contained in D and let 

U = {x e B r \V(x) > 0} (2) 

Additionally let V > throughout U. Then the point x = is unstable. 

Here, the instability referred to, is instability in the Lyapunov sense. When 
translated in terms of the particle picture, this amounts to exactly the same 
sense as conveyed in [3]. However, this result is more general and can be ap- 
plied even to a classical field theory, like General Relativity, which cannot be 
quantized unambiguously at present. Now, we apply this result to the case of 
the system with a Lagrangian that depends non-degenerately on the second and 
higher order derivatives of the dynamical quantities. For simplicity and without 
loss of generality, we treat the case of a system with a finite number of degrees 
of freedom. 

Ostrogradsky Instability and its consequences 

In [3], the Ostrogradsky Hamiltonian for one degree of freedom has been con- 
structed explicitly for a higher derivative Lagrangian L(q, q, q) to be of the 
following form 

H{qi,q2,pi,pi) =Pi<72 +KQi > 32,02) - L(q 1 ,q 2 ,n) (3) 
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Here, the choice of canonical coordinates is as follows [3 



d r d dL 



(4) 




(5) 



An approach to proving the instability of a system was already carried out 
in [7], in the context of charged solitons. The ideas presented below are a 
combination of the steps taken in that paper together with the result quoted 
above. Now, we come to the crux of the result. Phase space translations are 
canonical transformations [5]. Therefore, by a suitable canonical transformation, 
the Hamiltonian can be thrown into the following form: 



Here, h(Q, P) is a function that contains no linear terms. Therefore, in a neigh- 
bourhood of the origin, the first term dominates. If we choose the Chaetev's 
function as V(Q, P ) = Q1Q2 , we get for the derivative of this function along 



Again, the function u(Q, P) is dominated by the first term in a neighborhood of 
the origin. Therefore, in a small neighborhood around the origin, the conditions 
of Chaetev's theorem are satisfied so that this theory is unstable at the origin. 
Finally, by a suitable translation, we can reach the same conclusion at any point 
of the trajectory. Hence, we conclude that such a theory is unstable everywhere. 

Now that we have shown that any higher derivative Lagrangian is unsta- 
ble in the Lyapunov sense, we can return to the question of modifying gravity 
to explain accelerated expansion. In the literature, f(R) modifications are dis- 
cussed very extensively, where R is the Ricci scalar [3] , [9] . It is claimed that no 
other non-trivial modification can satisfy the constraints imposed by the result 
proved above [3]. However, it turns out (as shown below) that Ostrogradsky's 
theorem can be overuled for Born-Infeld (B-I) type of modifications to the E-H 
Lagrangians. 

Lagrangians are particularly attractive because of their intriguing properties 
[10] . The B-I Lagrangian first made its appearance with the modification of the 
electrodynamic action by Born and Infold. This was done to remove the infinite 
self energy of point charges by introducing an upper bound on the magnitude 
of the electric field. It was further pointed out that the B-I electromagnetic 
field propagates without birefringence [10]. Analogous work has been carried 
out for the case of the gravitational field. The relevant equations of motion can 
be obtained by varying the action given below : 



H(Q,P) = Q 1 P 2 + h(Q,P) 



(G) 
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Two dynamical quantities enter the action in (8). They are the metric and the 
connection. The Palatini formalism results if these quantities are considered 
independent degrees of freedom. The action is varied with respect to both the 
metric and the connection. The more usual method of varying gravitational 
actions explicitly uses the Levi-Civita connection and treats the metric as the 
sole dynamical freedom of the gravitational field. Here, we briefly recapitulate 
(after [5]), the equations of motion resulting from (8) 

= 9^ + kR^ (9) 



V\9 



q» v = \g» v - kT^ (10) 



Note the slight subtlety in notation. q^ v represents the inverse of q^. A is a 
cosmological constant which is necessary for consistency. The variation with 
respect to the connection yields 



pa 



The cosmology based on the Born-Infeld Theory has been restricted to the 
case where the spatial hyperspaces are flat [11]. We have extended these results 
to the case where the spatial hyperspaces can be curved. Before proceeding 
further, we must mention here an important consequence of the Ostrogradsky 
instability which shows why the Palatini formalism is important for the B-I 
lagrangian. 



The Born Infeld Lagrangian can only be varied in the Pala- 
tini formalism 

Suppose (8) is varied in the metric formalism. Then as stated in above, the 
connection is fixed as the Levi Civita connection. In that case, the determinant 
in (8), can be expanded in terms of the trace as follows (upto third order) 



det{g„ v + kR„ v ) m 1 + - {kR - kK - k 2 S) (12) 

where K = R^R^ - \R 2 and S = m^R^R* - 6RR^R^ + R 3 . 
Therefore, if the metric is considered the sole dynamical quantity, the La- 
grangian would involve the second and higher derivatives of the metric. But 
this would result in an unstable theory in accordance with the result proved 
above. 

However, since only the first derivatives of the connection coefficients appear 
in the expression for the Ricci tensor, there is no instability provided the ac- 
tion in (8) is varied by taking the metric and connection as separate degrees of 
freedom. Therefore, Born-Infeld Lagrangians can be varied consistently in the 
Palatini formalism. 
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The FRW equations in Born Infeld Gravity 



With the question of viability settled, we can now set up and solve the FRW 
equations for B-I gravity. We generalize here the flat space results of [11] to 
curved spatial hyper surfaces As usual, we assume that the universe is ho- 
mogenous and isotropic. This results in the following ansatz for the metric in 
comoving coordinates [2] 

9llv dx»dx v = -(dx ) 2 + a{tf | + r 2 {{d&) 2 + sm 2 0(#) 2 ) J (13) 

For the B-I equations, we also need an ansatz for the auxiliary metric q^ v . This 
is at once supplied by the demand for homogeneity and isotropy as [6] 

<Zoo = -U{tf (14) 
Qii = V(t) 2 %J (15) 

Here, 7^ represents the spatial components of the metric. Equations (14), (15) 
and (11) then determine the non-zero connection coefficients as 



U t = [ a 2 VV adV 2 \ 

U rr \ (1 - kr 2 )U 2 (1 - kr 2 )U 2 J 



a 2 VV adV 2 



U 2 U 2 



= TLsin 2 9 



d V 
a + V 



= -sinOcosO T^ g = cote 

rg r = rj r = i r^ = r^sin 2 

The [w — 00 and 11 components of the Ricci tensor are given by 
Roo = SH!(H + H 2 ) - 6HH 2 - 3 • 

+ V 2 ^H 2 + ^ + VV + HV 2 {H + H 1 -H 2 }}} 

Here, H = H 1 = %, H 2 = % 

For the matter distribution, we assume a perfect fluid with pressure p and 
density p. It can further be shown that the fluid satisfies the equation of con- 
tinuity. Finally, we can write the following equations of motion using (9) and 
(10) ' 



G 



1 - U 2 = k ^SH^H + H 2 ) - 6HH 2 - 3 
2kn 



a V 
a + V 



V 2 -l = ^ + -^\2HVV + V 2 

+ VV + V 2 |iJ 2 + ^ j +VV + HV 2 {H + Ht -# 2 }} 



UV = A + up 
V 3 



u 



A + up 

The fluid also has equations of state and continuity 

jp 

-3H(l + w) 



p = wp 

P 
P 



After doing all the manipulations, we arrive at our final expression for the 
Hubble rate in B-I cosmology 

. (^5 + {4K + 2/3 + 2 + a (i a} {(£)2 - 1) - 2k) + { -^}} a 

H ~a~ 2(a 2 + 2f3 + 2 + a@ - a) 

(16) 

Here rv — 3K P( 1 + t ") f m > 1 ] anr | « — 3«p(l+ m ) f 3m 1_ \ 

nere, a- 4 | A +Ktop + x+wp j ana p - 4 | a+ k ?«p A+^pJ- 



As illustrated in Figure 1 and 2, the main point of difference between the 
solutions of the FLRW equations in Einstein's gravity and in the present case is 
that closed universes in B-I gravity do not experience a 'big-crunch' even after 
the elapse of a large amount of time. Moreover, the radius always accelerates in 
this theory, which can be contrasted with the case of general relativity, where it 
always decelerates. An important point that was observed during the plotting 
of these figures was that the actual values of the radius as a function of time de- 
pended critically on the values of the coupling constant k and the cosmological 
constant A. Therefore, by tuning these values properly, one can, in principle, 
make the actual values of the radius agree with the experiment. 

On the other hand, an advantage of f(R) modifications is that any history of 
the evolution of the universe can be supported by a proper choice of the function 
f (R) [3] . This flexibility has been lost in the case of the Born- Infold modifica- 
tion. Nevertheless, as the plots show, the history predicted by our calculations 
is to a tolerable extent, consistent with observations. 
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X 10~29 seconds 



Figure 1: The Hubble rate as a function of time for a matter dominated open 
universe. Present time corresponds approximately to abscissa = 0.01 




X 1Q-29 seconds 



Figure 2: The Hubble rate as a function 
of time for a matter dominated closed uni- 
verse. Present time corresponds to abscissa 
= 0.01. A comparison with the previous fig- 
ure indicates that the Hubble rate increases 
much faster for a closed universe than an 
open universe. 




X 10-29 seconds 



Figure 3: The cosmic radius as a function 
of time for a matter dominated closed uni- 
verse. As discussed in the text, the acceler- 
ation of the universe is clearly apparent 
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